PSEUDO LINEAR ATTITUDE DETERMINATION OF SPINNING 

SPACECRAFT 


Richard R. Harman* and Itzhack Y. Bar-Itzhack* 

This paper presents the overall mathematical model and results from pseudo linear 
recursive estimators of attitude and rate for a spinning spacecraft. The measurements 
considered are vector measurements obtained by sun-sensors, fixed head star trackers, 
horizon sensors, and three axis magnetometers. Two filters are proposed for estimating 
the attitude as well as the angular rate vector. One filter, called the q-Filter, yields the 
attitude estimate as a quaternion estimate, and the other filter, called the D-Filter, yields 
the estimated direction cosine matrix. Because the spacecraft is gyro-less, Euler’s 
equation of angular motion of rigid bodies is used to enable the estimation of the angular 
velocity. A simpler Markov model is suggested as a replacement for Euler’s equation in 
the case where the vector measurements are obtained at high rates relative to the 
spacecraft angular rate. 

Extended Abstract 

q-Filter Dynamics 

The first dynamics equation we consider is the following Euler’s equation for the 
angular motion of a spacecraft (SC). It is [1, pp. 522, 523] 

ri> = I -1 [(Ico + h)x]oo + 1' 1 (T - h) (1) 

where I is the SC inertia matrix, e> is the angular velocity vector, h is the angular 
momentum of the momentum wheels, and T is the external torque acting on the SC. The 
symbol [ax] denotes the cross product matrix of the general vector a. Attitude is 
represented by the attitude quaternion whose kinematic equation is [1, p. 512] 
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In the q-filter we augment Eqs. (1) and (2) to form the following dynamics 
equation, which includes the noise terms 

ni r mi r t-> fr _ i\i Tw i 

(4.a) 

The unbiased white-noise vector w m accounts for the inaccuracies in the modeling of the 

SC angular dynamics, and w q is an unbiased white-noise vector that accounts for 

modeling errors in the quaternion kinematics. 

When the measurements come at a relatively high frequency we may be able to 
replace the* SC angular dynamics model in Eq. (4.a) with a simpler Markov model [2]. 
Consequently, Eq. (4.a) is replaced by the model 
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where [t] is a diagonal matrix whose elements are the inverse of suitable time constants. 


q-Filter Measurement Model 


b jm =[0 3 HjOj.q)] 
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where 
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rj is the reference vector corresponding to vector sensor j, and v jb is white noise. 


D-Filter Dynamics 


Using Euler’s equation and assuming the spacecraft attitude is represented as a direction 
cosine matrix, the dynamics take on the following form: 
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where d T = [df d 2 


d 3 ] , dj T is the transpose of the jth column of the direction cosine 


matrix, and CD = 


[d,x] 
[d 2 x] 
[d 3 x] 


where [dj x] is the skew symmetric matrix for jth column of the 


direction cosine matrix. 


The D-Filter Measurement Model 

For vector measurements, b jm = [d,r, | d 2 r 2 1 d 3 r 3 ] + v jb where r is the corresponding 

reference vector for the observation and dj is the ith column of the direction cosine 
matrix. This equation can be rearranged to form the measurement model: 
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Conclusion 


Both the q-Filter and the D-Filter will be tested against simulated data and a comparison 
will be made of the relative performance of each. 
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